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I. Introduction 

Engineering the artificial capture of asteroids in the neighbourhood of the Earth is a challenging 
problem. The simplest method would be to alter an asteroid’s Jacobi constant so that it is trapped 
within the zero velocity curves (ZVCs) in the neighbourhood of the Earth. This has been investigated 
by Q who considered a range of NEOs in the context of the Earth-Sun Circular Restricted Three- 
Body Problem (CRTBP) and determined the Av needed to close the zero velocity curves and truly 
capture the object. The lowest Av required was found to be of the order of 400 m s -1 . Another 
method suggested by 2j uses a variant of patched conics to transfer an asteroid into the Earth’s 
Hill sphere combined with an impulsive transfer onto a bound geocentric orbit. They identify 23 
asteroids that could be captured within 10 years, which have final capture Av requirements ranging 
from 700 m s -1 to 4900 m s -1 . Similarly, 3] and 4] consider various capture strategies involving 
patched conics to create resource maps of accessible NEOs within a given Av range. 

An alternate method for capturing asteroids has been suggested by |5j] who find a new class of 
‘Easily Retrivable Objects’ (EROs) among the near-Earth object (NEO) population. These objects 
are those that can be easily transferred onto the invariant manifolds associated with Lyapunov 
and Halo orbits around the Earth-Sun L\ and L 2 points. They find a list of 12 objects that can be 
transferred into the Earth locale with a An of less than 500 m s^ 1 . The problem of finding long term 
stable parking orbits and associated transfers for such asteroids has received less attention though. 
In 6], families of linearly stable co-orbital periodic orbits about a central body are identified. Such 
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orbits are potential locations for ‘parking’ asteroids about an object such as the Moon, although 
transfer strategies to such orbits are not immediately apparent. 

A different approach to the problem is used by 7], who consider ‘temporarily captured asteroids’ 
(TCAs) such as 2006 RHi 2 o- TCAs are objects that temporarily come within the Hill sphere of 
the Earth, and remain for a short period of time (months to years). Using a dynamical model of 
the Earth-Moon-Sun system they find that a relatively low Av of 32 m s _1 applied gradually over 
several months can extend the capture phase by several years. 

Another method is introduced by jj|, who considers the result of close approaches of binary 
asteroids with the Earth, demonstrating that under some conditions a binary exchange can occur. 
The result of such an encounter is the capture of one of the asteroids about the Earth, while the 
other is ejected from the system on an escape trajectory. 

In developing these new strategies for asteroid capture, astronomical theories explaining the 
capture of irregular satellites of the giant planets are relevant. There are three main theories: gas 
drag (capture due to dissipation induced by the gas drag in the early solar system), pull down cap¬ 
ture (capture as a protoplanet’s Hill sphere expands as its mass increases) and the broad category 
of n-body interactions. Of the first two mechanisms pull down capture does not readily translate 
to methods for engineering asteroid capture, while ideas relating to gas drag have been considered 
by |9( who propose a method of aero-assisted capture that exploits the dissipation of energy due to 
atmospheric drag to capture asteroids Earth-bound orbits. The last category presents further oppor¬ 
tunities, and has been the focus of several studies already, such as the binary exchange mechanism 
proposed by 8] or the invariant manifold method of 5] both discussed above. 

A review of the various n-body interaction theories in the context of formation of irregular 
satellites is presented in m These fall into the broad categories of collisions, n-body interactions 
due to resonance crossings of the giant planets and exchange reactions, such as the disruption of a 
binary asteroid. The disruption of a binary asteroid has been investigated in the context of asteroid 
capture by 8], but presents engineering challenges as it is very sensitive to the configuration of the 
asteroid during the encounter. 

Another dynamical mechanism for the capture of irregular satellites is chaos-assisted capture, 


2 
















as proposed by 111. Here the phase space of the Sun-Jupiter CRTBP is shown to consist of islands 
of Kolmogorov-Arnold-Moser (KAM) orbits separated from scattering orbits by a thin chaotic layer. 
A small perturbation such as dissipative gas drag can then capture a chaotic object into this regular 
region of phase space. In a two-dimensional model the captured object can then not escape, while 
in three dimensions Arnold diffusion means it can move out through the KAM tori exponentially 
slowly, and is thus effectively captured on astronomical timescales. They show that this mechanism 
is viable in the context of both the two-dimensional and three-dimensional Sun-Jupiter CRTBP. 
In 12j, the chaos-assisted capture mechanism in the Eccentric Restricted Three-Body Problem 


(ERTBP) is shown to be robust to perturbations from other planets as well. In 


13| the formation 


of Kuiper-belt binaries using this capture mechanism is studied in the context of the three- and 
four-body Hill problem. In this case a temporary chaotic binary is moved to a regular trajectory 
through scattering encounters with other objects. The chaos-assisted capture mechanism provides 
an interesting possibility for engineering a relatively low energy artificial capture of an asteroid in 
the Earth’s neighbourhood. 

In this work we show how a temporary captured asteroid that passes close to a regular region 
of phase space can be moved onto the nearby KAM tori and essentially be captured permanently 


in the Earth’s neighbourhood, via the chaos-assisted capture mechanism of 


ll| . We demonstrate 


the method and the relatively low Av required for an example asteroid trajectory. 

We use the two-dimensional Hill problem as a simple first model of the Earth-Sun system, and 
calculate trajectories and Poincare surfaces of section (SOS) using a 6th order Runge-Kutta-Nystrom 
(RKN) symplectic integration scheme. Although this does not capture the three-dimensional dy¬ 
namics or the effect of the Moon on the near-Earth environment, it does demonstrate the basic 
principles of the method and provides a demonstration of the capture energies required. Further 
development is not expected to vastly alter either of these, merely the details of the trajectories. 
Details of the mathematical model and integration scheme are given in the next section, followed 
by the results for the example trajectories. 
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II. Hill problem surface of section 


The Earth-Sun system can be modelled simply using the two-dimensional Hill problem. To 
generate a Poincare surface of section (SOS) numerical methods are needed to integrate the equations 
of motion. One such method is a symplectic integration scheme. Such integrators preserve the phase 
space structure and so are advantageous for the generation of Poincare SOSs. A 6th-order symplectic 
integration scheme for the Hill problem can be easily constructed by extending the method of 14|, 


who derive a second-order symplectic scheme for the planar Hill problem, to a higher order scheme 


such as given by [15]. The planar Hill problem has Hamiltonian 


1 


H{x, y,p x ,p y ) = P 2 x +P 2 y ) + {yp x ~xp y )-x 2 + ^y 2 - - + /2 


(1) 


where x and y are the scaled coordinates relative to the Earth in the rotating frame and p x = x — y 
and p y = y + x are the corresponding conjugate momenta. This is the dimensionless formulation of 
the Hill problem and uses units such that the mean motion of the Earth about the Sun is unity. It 
can be derived from the CRTBP by moving the origin to the location of the Earth and taking the 
limit of the mass ratio going to zero. In the Hill problem the energy E is a constant of the motion 
and is often used in the form of the Hill-Jacobi constant Cr = —2 E. In 14], a canonical transform 


is introduced that simplifies the Hamiltonian and allows it to be split into two analytically integrable 
parts. The new coordinates and momenta are given by X = x, Y = y, Px = p x + y , Py — Py + x 
and the transformed Hamiltonian is 


H(X, Y, P X ,P Y ) = \(P‘x+ Py) - 2 XPy + ^ - (;g2+ y 2)1/2 ■ 


( 2 ) 


This Hamiltonian can be split as H = Hi + H 2 with 


Hi = \{P%+P$)-2XPy 

(3) 

X 2 ^ 

H2 2 (X 2 + Y 2 ) 1 / 2 ' 

(4) 
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The equations of motion associated with both H\ and Hi are now trivially integrable. As given by 


14j the evolution from time t to t + r under Hi is given by 


X(t + r) = X(t) + rPxit) + T 2 P Y {t) 
Y(t + t) 

Px(t + r) 

Pv{t + T) 


Y(t) + r(Py(t) - 2 X(t)) - r 2 P x (t ) - -T 3 Py(t) 
P X (t) + 2rPy(t) 

Py(t) 


(5) 


and the evolution from time t to t + r under Hi by 

X{t + T) = X(t) 

Y(t + r) = Y(t ) 

Px(t + T) = Px{t) - T (X(t) + 
Py(t + T) = Py{t)-T 


m 


(. X(t ) 2 + Y{t) 2 f/ 2 
Y(t ) 


(X(t,) 2 +Y(t,) 2 ) 3 / 2 J ' 


( 6 ) 


In 


14l |. the authors further partition the term Hi into two pieces to permit quick detection of 


collisions between particles for simulations of ring dynamics. However, this is not required here 
and the splitting above is sufficient. Further, this splitting H = Hi(P,X) + Hi{X) satisfies the 
condition to use a Run ge-K utta-Nystrom (RKN) symplectic method, as the term Hi is quadratic in 
the momenta (see e.g. 


(also used by 


15J). Such a method is the 6th order RKN symplectic scheme given by 


15] 


16( for the Hill problem in regularized coordinates). We use a 6th order method for 


accuracy as the two Hamiltonians Hi and Hi are fairly comparable in size. The RKN symplectic 
scheme is more accurate than other schemes at this order, but requires more force evaluations per 
step. However this is not an issue here as we use a symplectic algorithm for accuracy as explained 
above, rather than speed. The RKN scheme was found to be sufficient to generate the SOS needed 
for this work. 

Let (j>Hi represent the mapping associated with the Hamiltonian Hi and (pni that with Hi. The 
RKN-6 scheme approximating the evolution under the Hamiltonian H = Hi + Hi over timestep r 
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is 


<I>rkn6(t) = </>#! (aiT)<j>u 2 ( biT)(j> Hl (a 2 T)4>H 2 {b 2 T)<f> Hl (a 3 r)</) ff2 (b 3 T)4> Hl (a 4 r)0 ff2 (6 4 t) 

<^H! (a 4 r)^) ff2 (b 3 r)(j) Hl ( a 3 T)4>H 2 ( b 2 r)4> Hl ( a 2 r)(f> H2 (bir)(j) Hl (air)^ (6ir) (7) 

where the coefficients a,- and 6, are 


ai = -1.01308797891717472981 fo 4 = 0.00016600692650009894 


a 2 = 1.18742957373254270702 & 2 = -0.37962421426377360608 


a 3 = -0.01833585209646059034 6 3 = 0.68913741185181063674 


a 4 = 0.34399425728109261313 6 4 = 0.38064159097092574080 


( 8 ) 


as given by 


13- 


The SOS is calculated as follows. The SOS chosen is the plane in x-y defined by p x = 0 and 
with y > 0 (see the following Section for more details on this choice) i.e. it is a plot of the points at 
which trajectories cross the x — y plane with zero a:-momentum and positive y-velocity. 

A trajectory is generated by selecting an initial point on the section and using the constraints 
p x = 0 and y > 0 along with the energy integral to calculate the initial momenta of the orbit. The 
initial state is then propagated forwards in time using the symplectic integration scheme described 
above. Intersections of the trajectory with the SOS are checked for after each timestep by looking 
for a change in sign of p x to detect the crossing of p x = 0. If a crossing is detected the exact 
location of the intersection is then interpolated. With a small enough timestep this provides an 
accurate calculation of the SOS. The integration length was chosen to be sufficient to fully populate 
the KAM tori. For the results presented here a timestep of 0.005 in dimensionless units is used 
(approximately 6 hours) and the integration length is approximately 0.8 MYr. The method was 
tested by computing specific sections that could be compared to results in the literature. 

The initial points for each SOS are selected using an interactive GUI that generates the trajec¬ 
tories from points selected with the mouse. This allows individual features on the surface to be fully 
investigated, which could otherwise be missed if a random set of initial points was used instead. 
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III. Chaos-assisted artificial capture of an asteroid 


The SOSs of the planar Hill problem are well studied, see for example 17] and 13]. As found 
in these studies the phase space near the central object (the Earth here) at energies at above that 
of the two Lagrangian points consists of KAM regions surrounded by thin chaotic layers embedded 
in regions of escape trajectories. We select at random an example initial chaotic asteroid trajectory 
that enters the Hill sphere, passes near to a region of KAM tori surrounding a prograde periodic 
orbit (one of the g' family) and subsequently escapes the Earth neighbourhood. The trajectory is 
somewhat arbitrary, but sufficient to demonstrate the method and the energy requirements. The 
energy level of the initial asteroid orbit is E = —2.152 (in Hill units), which is similar to the examples 


given by 13] for the Hill problem for binary asteroids. At this energy level the ZVCs are open at 
the Lagrangian points, but the Hill sphere is still a distinct region separated from the rest of the 
phase space by excluded areas. 


The example initial trajectory is shown in configuration space in Figure 1(a) with the corre¬ 
sponding surface of section (SOS) in Figure [5] The SOS used is the plane x-y that has p x = 0 and 
y > 0, as used by ll| ]l3] for the planar CRTBP and planar Hill problem. The SOS clearly shows 
the KAM region surrounded by a thin region of chaotic motion. Trajectories that pass into this 
region of chaos are sometimes said to be ‘sticky’, as they stay within it for a number of crossings 
of the SOS, often ‘sticking’ to tori embedded within it. The behaviour can be seen in the x-y plot 


of the trajectory in Figure 1(a) the object enters the Hill sphere, roughly follows the shape of a g 1 
family orbit for a number of revolutions, before exi ting the sphere once more. 


The natural chaos-assisted capture method of 


11] theorises that orbits similar to this can be 


moved into the nearby regular region of phase space (at a different energy level) by effects such 
as gas drag. We propose a method for capturing asteroids that does this through an artificial 
instantaneous Av applied to the asteroid. Essentially, the chaotic orbit will be ‘kicked’ to a new 
energy level where the trajectory is now regular. Using the SOS for nearby energy levels we can 
identify points at which the original trajectory crosses the section nearby to a regular region. The 
SOSs of nearby energy levels are shown in Figure [3] As noted by 11], depending on location either 
an increase or decrease in energy may move the asteroid into a regular KAM region of phase space, 
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and as such the SOS has been plotted for —2.154 < E < —2.150. It can be seen that the regular 
region changes fairly significantly with even a small change in energy, and some of the crossings of 
the SOS of the original asteroid trajectory which are chaotic at E = —2.152 are within a regular 
region of phase space at E = —2.1535 and E = —2.154, as well as at E = —2.1505 and E = —2.15. 
An instantaneous impulse at one of these points is all that is required to transfer the chaotic orbit 
to a regular orbit. 

Taking the impulse to be perpendicular to the SOS at the moment of transfer allows the Av to 
be calculated for each potential transfer point. The transfer at E = —2.1535 is the best, requiring 
Av = 64.9 m s -1 in the negative y direction. The transfer at E = —2.1505 is slightly less favourable, 
requiring Av = 105.5 m s -1 in the positive y direction this time. By comparison, to close the ZVC 
at the first point through a similar impulse would require a Av of approximately 500 m s -1 (and 
approximately 840 m s -1 at the second point). It is clear that the new method leads to significant 
savings in capture Av cost. 

The exact point of the best transfer is marked on Figure [3(b)| and the new trajectory shown in 


Figure 1(b) With a relatively low energy requirement the asteroid has been essentially permanently 
captured in a prograde orbit about the Earth. It has the further advantage of being protected from 
collision with the Earth as well as escape from the Hill sphere. 

The method used here is simple and merely intended to demonstrate the concept. To extend 
the work to more realistic models a three-dimensional dynamical model of the Earth-Moon-Sun 
system is essential. An optimization scheme combined with a method to more precisely determine 
the re gion s of chaotic and regular trajectories, such as the fast Lyaponuv indicator (FLI) maps used 
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13| would provide a more robust means of detecting the optimal transfer of an asteroid using 


the chaos-assisted capture method. However, the results presented here demonstrate that even using 
a simple SOS method to select the transfer point yields very favourable Av requirements. 

Finally, the transfer scheme is limited in that it requires the asteroid to be in the chaotic layer, 
however we note that other authors, such as [7j have looked at transferring NEOs onto these types of 
trajectories. The artificial chaos-assisted capture method can therefore be seen as the final capture 


stage in a more general strategy for transferring asteroids into the Earth neighbourhood. 
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(a) The initial trajectory and ZVC at E = —2.1520 (b) The new trajectory at E = —2.1535 

Fig. 1 The initial and final example trajectories in configuration space. The corresponding 
ZVC surface is also shown in gray in each case (the ZVC closes at E ~ 2.1635) 
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Fig. 2 The surface of section for the initial trajectories energy level E = —2.152. The crossings 
of the example trajectory are plotted in red 

IV. Conclusion 

The theory of chaos-assisted capture of irregular satellites of the giant planets can be used to 
develop low-energy capture strategies for certain classes of Near Earth Objects (NEOs) that pass 
close to a chaotic layer associated with a regular Kolmogorov-Arnold-Moser (KAM) region of phase 
space. Investigation of the surfaces of section (SOS) within the two-dimensional Hill problem model 
of the dynamics can be used to identify points along a trajectory at which a low energy instantaneous 
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(f) E = -2.1510 
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(g) E = -2.1505 (h) E = -2.1500 

Fig. 3 The surface of section for the a range of energies close to the original energy level 
E = —2.152. The top four figures are at a lower energy and the bottom four figures at a higher 
energy. The locations of the SOS crossings at the original energy of the initial trajectory are 
shown in red on each plot. The best potential transfer point is identified as a black point on 

the E = —2.1535 plot 


transfer can be made to move the object into the regular region and onto a long-lived and essentially 
bound orbit about the Earth. The strategy places the captured object in a regular region of phase 
space that, in the simple model used, is protected from both collision with the Earth and escape 
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from the Hill sphere. The method of transfer does not require a control strategy, and further neither 
does the final orbit require any control mechanism in this model. The Av requirements are low, 
comparable to those needed for temporary capture strategies and much lower than those needed by 
other permanent capture strategies. Here a Av of order of tens of m s^ 1 is required, compared to 
the hundreds to thousands of m s _1 needed for other capture strategies. 

The method is limited in that it is only applicable to asteroids that pass close to a chaotic layer 
near a regular region of phase space. However, transferring an asteroid artificially into this region of 
phase space has been addressed by other authors and studies, who consider low energy strategies for 
transferring asteroids onto longer lived (but still temporary) ‘sticky’ chaotic orbits. Given this, the 
chaos-assisted capture scheme presented here provides a method for the final stage in the artificial 
capture of an asteroid. 

We also note that the use of an instantaneous impulse to transfer the asteroid to the regular 
region is likely to be unworkable for any asteroid of significant mass, it is merely a convenient 
method to illustrate the magnitude of the change in energy required to capture an asteroid using 
this method. A sustained period of low-thrust impulse would be a more practical solution, and is 
in fact more in analogy with the original theory of formation of irregular satellites, which suggests 
that they are transferred into the regular KAM regions via weak dissipation occurring over a period 


of time 


111 . 


The example given in this work is a relatively low energy orbit and here the SOS changes fairly 
rapidly with energy. For high energy objects the SOS changes slower so the energy requirement 
and Av will necessarily be larger. However, it worth noting that the true capture Av for such 
objects can be in the region of tens of thousands of m s -1 so the saving in energy is still likely to be 
significant. It is well known that at higher energies the regular region associated with the prograde 
orbits no longer exist and the transfer has to be to regions surrounding retrograde orbits. The phase 
space here contains trajectories that come close to collision with the primary and we note that this 
provides an interesting opportunity to deflect and capture objects that would otherwise collide with 
the Earth. 

This work is intended to illustrate the potential chaos-assisted capture as a strategy for the 
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artificial capture of asteroids. Extension to a more realistic three dimensional model that includes 


important perturbations such as the effect of the Moon, as well as full optimization to select the 
lowest possible energy transfer is necessary. However, even without this the results demonstrate 
that this is a viable method for the final capture stage of asteroids into the neighbourhood of the 
Earth. 
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